Interplay between bending and stretching in carbon nanoribbons 
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We investigate the bending properties of carbon nanoribbons by combining continuum elasticity 
theory and tight-binding atomistic simulations. First, we develop a complete analysis of a given 
bended configuration through continuum mechanics. Then, we provide by tight-binding calculations 
the value of the bending rigidity in good agreement with recent literature. We discuss the emergence 
of a stretching field induced by the full atomic-scale relaxation of the nanoribbon architecture. We 
further prove that such an in-plane strain field can be decomposed in a first contribution due to the 
actual bending of the sheet and a second one due to the edges effects induced by the finite size of 
the nanoribbon. 
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PACS numbers: 62.25.-g, 62.20.D-, 46.70.Hg 



I. INTRODUCTION 

Graphenei plays a unique role in materials science 
since it is the mother structure of most carbon sp^ 
nanosystems of current interest. By stacking, folding 
or bending a graphene sheet it is indeed possible to 
generate, respectively, graphite-like systems, fullercnc 
cages (pentagonal rings arc here needed as well) or 
nanotubes. In particular, the bending properties are 
critical in attaining the structural stability and mor- 
phology for both suspended and supported graphene 
sheets, and directly affect their electronic propertiesi^ 
Rippling of pure graphene has been also observed with 
mcsoscopic amplitude and wavelength, both for sus- 
pended monolayers^ and sheets deposited on substrates 
such as silicon dioxide^ Moreover, the bending prop- 
erties play a central role in the design of graphene- 
or carbon nanotube-based devices, like e.g. mechani- 
cal resonatorsi^ The bending features of functionalized 
graphene sheets have been probed by atomic force mi- 
croscopy, observing that the folding behavior is domi- 
nated by defects and functional groupsj^ Finally, bend- 
ing ultimately governs the carbon nanotubes unzipping 
process, recently used to produce narrow ribbons for 
nanoelectronics.— With the same technique, a new class 
of carbon-based nanostructures, which combine nanorib- 
bons and nanotubes, has been introduced in order to ob- 
tain magnetoresistive devices 

Within this scenario we frame the present investiga- 
tion, addressed to improve our fundamental understand- 
ing of the bending properties of a one-atom thick carbon 
sheet. The main goal is twofold: (i) to draw a thor- 
ough theoretical picture on bending, fully exploiting the 
elasticity theory and providing an atomistic quantitative 
estimation of the corresponding bending rigidity; (ii) to 
prove that the bending process of a carbon nanoribbon is 
always associated with the emergence of a (small) stretch- 
ing, particularly close to the edges. These results have 
been obtained by combining continuum elasticity theory 
and tight-binding atomistic simulations (TB-AS). 

The conceptual development and actual exploitation 



of our theoretical model proceeds through the following 
steps. At first, by means of continuum mechanics we have 
obtained the exact shape for a purely bended nanoribbon, 
by imposing suitable boundary conditions. The bending 
rigidity is then evaluated by TB-AS for several nanorib- 
bons differing by length and width. As a second step, we 
observed that, under the above assumption of pure bend- 
ing, the corresponding rigidity must be a constant inde- 
pendent of the actual shape of the sheet. Nevertheless 
by allowing full atomic-scale relaxation during bending, 
we rather found a geometry-dependent rigidity, a feature 
that we have attributed to the onset of stretching phe- 
nomena. Therefore, as final step, we have developed a 
procedure to discriminate between stretching and bend- 
ing energy, so providing a complete picture about the 
mechanical behavior of graphene and also reconciling the 
atomistic data with the continuum theory results. 

The structure of the paper follows: in Section II we 
outline the theoretical framework from both the con- 
tinuum elasticity theory and the tight-binding atomistic 
simulations point of view. In Section III we describe the 
results concerning the bending stiffness and the interplay 
between stretching and bending. Finally, in Section IV 
we draw the conclusions. 



II. THEORETICAL FRAMEWORK 

A. Continuum picture 

The graphene strain energy density U [cVA~^] is de- 
fined asi^ii^ 



(1) 



where E [Nm^^] and v are the two dimensional Young 
modulus and the Poisson ratio, while k [eV] and R [eV] 
are the bending rigidity and the Gaussian rigidity, re- 
spectively. The in-plane deformation (stretching) energy 
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[given by the first two terms in Eq.(IT|)] is described by 
the standard small strain tensor e = i(VM + Vu'^), be- 
ing u the displacement field. On the other hand, the 
out of plane deformation (bending) energy [given by the 
last two terms in Eq.((T])] is described by the mean cur- 
vature % = -^li^ [m~"^] and by the Gaussian curvature 
K ~ kik2 [m~^], where ki and k2 are the principal curva- 
tures at a given point on the surfaceji^ They are straight- 
forwardly given by fci = I/-R1 and fc2 ^ I/R2 where Ri 
and i?2 are the principal radii of curvature at that point. 
In the case of a continuum plate of thickness h made 
of an isotropic and homogeneous material, the classical 
Kirchhoff theory provides k = j2Y^ ^^"^ ^ ~ I^rTi^ 
(note that E = Yh where Y is the three-dimensional 
Young modulus) For an infinitesimally thin graphene 
monolayer such a theory does not apply, since the thick- 
ness h cannot be unambiguously defined and the bend- 
ing moment has simply a different physical origin. While 
the bending moment for the Kirchhoff plate derives from 
a compression/extension of the different material layers 
forming the thickness h, in graphene it is due to the in- 
teractions among orbitals Pz which are affected by the 
bending process. Therefore, the determination of n and 
R for graphene is a well-posed (and, to a large extent, still 
open) problem, which is independent of the evaluation of 
E and 

Our model system is a rectangular ribbon with length 
I and width L (sec FiglT|) . The ribbon is bended without 
stretching {§ = 0) along its width. The boundary con- 
ditions consist in fixing the positions of the two parallel 
edges (with length ^) at a given distance a, while the at- 
tack angles is free to relax. This configuration involves 
only one curvature fci, leading to V. = ^ and K. = 0. 
By considering different values of a in the range (0, L) 
we obtained a set of differently bended configurations. 
The elastic problem consists in finding the sheet shape 
by minimizing the bending energy 



UdA ^\kI 



k\ds 



(2) 



where A = LI \s the total area of the system. If the 
configuration is described by the function z = z[x), then 



we have fci = z" / 1 + (z') 



3/2 



inhere z' = dz/dx and 



(P'zjdx^. On the other hand, ds = ^/gdx where 



^=^1 + [z'y. Therefore, Eq.® assumes the explicit 
form 



{z"f 



i + {z'r 



5/2 



dx 



(3) 



The problem consists in finding the curve z = z{x) min- 
imizing the energy functional in Eq.([3]) under the con- 
straint 



^Jl + [z'fdx 



= L 



(4) 




Figure 1: (Color online) Bended ribbon with length I and 
width L (red dashed line). The parallel edges with length I 
are fixed at distance a, while the attack angles 9 is free to 
relax. 



enforcing the absence of any in-plane stretching. 

By the application of the constrained variational calcu- 
lus we eventually obtain the final geometry in parametric 
representation [x{s),z{s)] 



x 
L 
z 



E((7)-f (am{K(g)(l-2^)},q 
cn{K(,)(l-2l)} 



K(9) 



(6) 



where s is the arc length {0 < s < L), q ~ sin | is the 
elliptic modulus and 6 is the attack angle given by 



a_oE(g) 



(7) 



The quantities E{q) and K{q) are the complete elliptic 
integrals, defined asi^J^ 



E(g)=J-(|,q), K(9)=£(|,q 



(8) 



where the functions T{v,q) and £{v,q) are incom- 
plete elliptic integrals of the first and second kind, 
respectivel y "'^^1-'^^ 



Tiv.q) 



S{v,q) 



da 



\/\ ~ q^ sin^ a 



1 — sin^ ada. 



(9) 



Moreover, by considering u ~ J-(v,q) we define the in- 
verse relation (with fixed modulus q) v = am{M}, which 
is called Jacobi amplitude function. Further, cn {u} = 
cosw = cos (am {m}) and sn{u} = sinf = sin (am {u}) 
are the Jacobi elliptic functionsJ^ Interesting enough, 
one can prove that liuia/L^o ^ = 130.709°, an universal 
value of the attack angle found whenever a = or L is 
very large. 
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Figure 2: (Color online) Bending rigidity k obtained for 
purely ribbons with several widths L. The average value is 
given by Kave = 1.4025 ± 0.0025 eV (yellow area shows the 
error bar) 



B. Atomistic simulations 

The present TB-ASiS have been performed making use 
of the sp'^, orthogonal, and next-neighbors tight-binding 
representation by Xu et aljii The present TB total energy 
model has been implemented within the scheme given by 
Goodwin et alr^ for the dependence of the TB hopping 
integrals and the pairwise potential on the interatomic 
separation. 

Applications to molecular-dynamics studies of liquid 
carbon and small carbon clusters indicate that this model 
correctly describes carbon systems over a wide range 
of environments 1^ This approach has been successfully 
used to show that the surface of nanodiamond particles 
reconstructs in a fuUerenelike manner, generating carbon 
clusters called bucky diamonds j2i Moreover, the growth 
of nanostructures (linear, ring, and fuUerenelike objects) 
in a carbon plasma^i and the formation of carbon clusters 
(onion-like and endohedral structures) from the conden- 
sation of liquid droplets^^ have been simulated by the 
present tight-binding model. Finally, this parametriza- 
tion has been recently used for determining the non- 
linear elastic moduli governing the graphene stretching 
elasticityJ^ 

The previous continuum analysis is useful both to cre- 
ate the input configurations for atomistic calculations 
and to define the simulation protocol. The investigated 
system consists in a nanoribbon formed by a perfect 
hexagonal carbon lattice, having width L in the range 
4-12 nm and length I imposed to obtain a simulation 
box containing a constant number of ~ 600 carbon 
atoms. Morever, periodic boundary conditions are as- 
sumed along the direction of the length /. The length 
(width) is developed along the armchair (zig-zag) direc- 
tion of the honeycomb lattice. Each nanoribbon is de- 
formed as defined in Eqs.([5]) and ([6|) in ten configurations 
corresponding to different values of a. In any bended 
configuration, all the interatomic distances are fixed at 
the equilibrium value for fiat graphene (so that no bond 
stretching is for the moment allowed) . The bending rigid- 
ity has been straightforwardly obtained as At = -^Ub with 
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Figure 3: (Color online) The theoretical results for the val- 
ues of the attack angles 8 (circles), predicted from Eq.(I7| for 
several ribbon with different width L, each at different edge- 
distances a, are compared with the corresponding data from 
atomistic simulations (crosses). 



Ub given by Eq.(l2]), where the integral 2 = kfcls is 
computed for the given configuration. It is important to 
remark that the obtained value for k must be indepen- 
dent of the actual configuration since the deformation is 
a pure bending one. 



III. RESULTS 

Accordingly to the scheme outlined in the previous Sec- 
tion, we have firstly evaluated the (pure) bending energy 
as Ub = £;;;™<*<=^ - i;/'"^*, where El"'^'^'"^ and repre- 
sent the TB-AS total energy of the bended (but not re- 
laxed) and equilibrium (flat) conflgurations, respectively. 
The atomistic results for k are reported in Figj^] (sym- 
bols) as function of the a/L ratio and for different width 
L. We estimate an average value Have = 1-40 eV. While 
the reported values of k (for nanotubes) vary in the range 
1 eV< K <2 eV^S we remark the most reliable ab-initio 
data K = 1.40 eVj^ and k = 1.46 eV^^ are in excellent 
agreement with our prediction, a feature standing for the 
reliability of the present computational procedure. 

Although reassuring, the above picture must be refined 
in order to properly take into account atomic-scale fea- 
tures. Therefore, full relaxation of the internal degrees 
of freedom of the bended systems is performed by zero 
temperature damped dynamics until interatomic forces 
resulted not larger than ~ 10~^eV/A. We have so gen- 
erated a new set of configurations where bending and 
stretching features are entangled. During the relaxation, 
the positions of the atoms belonging to the edges (i.e. 
atoms with a; = or a; = a, see Fig[T]) are fixed and. there- 
fore, the distance a between the edges remained constant. 
Overall we observed that the geometry is only marginally 
affected by relaxation as shown in Fig (31 Here we com- 
pare the attack angle 9 predicted from Eq.([7|) versus the 
ratio a/L with the corresponding values obtained from 
the relaxed configurations. We note that, for a/L — > 0, 
we obtain the universal value 130.709° as previously dis- 



4 




Figure 4: (Color online) Bending rigidity k computed by 
means of Ui = E'^^'^'^^itd- EHf^.^^^ (full circles). Straight lines 
correspond to the average value Kave = 1-40 eV as deduced 
from Fig.l. 



cussed. As a matter of fact, after the relaxation, the 
attack angle 9 do not change and the maximum varia- 
tion of L was as little as 0.005 nm, corresponding to a 
variation of the integral X smaller than 0.01%. Never- 
theless, even for such minor relaxations the energetics of 
the fully relaxed systems is expected to sizcably differ 
from the purely bended case, because of the extraordi- 
nary large value of the graphene Young modulus*^ It is 
therefore important to provide a new estimation of the 
bending energy for the fully relaxed configurations. 

Following the above argument, we evaluated the new 
bending rigidity n by means of the energy Uh = E'^tTaxed~ 
ElttLed and Eq.dU, where is the energy of a re- 

laxed bended ribbon and E^^l^^^ is the energy of a fiat 
ribbon after a full relaxation (different from the energy of 
the infinite graphene sheet because of the edge effects). 
In this case, we have found a variation of k upon a/ L 
as shown in Fig|3](full circles). This result suggests that 
atomic-scale relaxations upon bending have induced as 
expected an additional field of in-plane stretching, which 
provides new energy contributions as reported in Eq.(IT]). 
It is interesting to observe that the largest differences be- 
tween the unrelaxed and relaxed configurations are found 
for a/L ~ 1. In fact, in this case the forces exerted by 
the constraints (maintaining the distance a between the 
edges) are almost parallel to the graphene sheet, favoring 
the stretching emergence. 

This intriguing result opens the problem of how to 
disentangle bending and stretching features. As shown 
in FigUl this is especially important in the limit of 
small deformations, a situation of considerable practi- 
cal interest. To this aim we have defined a proof-of- 
concept computational procedure based on the virtual 
process of straightening (or unbending) of a given re- 
laxed and bended ribbon: atoms are projected from such 
a configuration onto a plane by conserving all the Inn 
bond lengths and all the 2nn planar angles. The pro- 
cess recovers a planar configuration, still maintaining 




Figure 5: (Color online) Comparison between the bending 
rigidity K computed through Ub = eI""'^'"^ - EI^"' (open cir- 
cles) and Ub = Ell^:^^ta - Eil:l,M.uea (crosses). The maxi- 
mum deviation is less than the 1.5%. 

all the details about any possible stretching (in-plane 
strain field); the corresponding energy £^ft"aig/jtene<i is 
straightforwardly evaluated by means of TB-AS. The 
bending rigidity n can be consequently determined by us- 
ing Ub = £^^eratfd - E^sirLghtened- rCSuhs are showu 

in FiglS] (crosses) where we also report k as obtained by 
Ub = E^^'^"'^'^'^ — eI'-"-* (open circles). The comparison 
points out a good agreement between the two different 
approaches since stretching features are either at all non 
considered (open circles) or included in both the bended 
and fiat configurations (crosses) so as to compensate. It 
is interesting to note that the constant trend of k versus a 
and L has been found similar to Fig|5] In other words, we 
have proved that the evaluation of k through the energy 

term Ub = KlTaxed' ^IttLed is not correct since it is cor- 
rupted by a strain energy amount which is not directly 
related to the bending process. The energy due to the 
sole stretching field (induced by the bending process) can 

be accordingly defined as E^f^^g^^^^^^ - i^/^^lerf- The 
demonstration that such an energetic contribution corre- 
sponds only to stretching relies on the fact that both the 
terms E^^^l^g^^^^^^ and ^^/^^^l^d have been evaluated on 
flat ribbons through TB atomistic simulations. 

A further evidence of the stretching emergence can be 
derived from FigIS] where the strain is calculated along 
the arc of length L (corresponding to the dashed line 
in FiglT]), labeled by the coordinate s. We can calculate 

three strain fields ei''J , si's* and ei*i which are respectively 
defined as the relative difference between: (i) the relaxed 
and straightened configuration (energy eH^^.^^^^^^^) and 

the flat relaxed configuration (energy -E'/ez'axed)' (ii) 
fiat relaxed configuration (energy E^l^^^^^^) and the flat 
unrelaxed configuration (energy £'/'"*); (iii) the relaxed 
and straightened configuration (energy E^^^^-gf^^^^^^) and 
the flat unrelaxed configuration (energy £'/'"*). While 
the strain eis^ is only due to bending, the term eit^ is 
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Figure 6: (Color online) eis' (strain induced by the bending), 
Ess' (strain induced by the edges) and ei'^ (total strain) versus 
s (red curves) . The gray scale map in background represents 
the same quantities in the sy-space for L = 12 nm and a/L — 
0.95 

induced by the presence of the edges (finite nanoribbon) 

in a flat configuration. The quantity represents the 
total strain induced by the relaxation of the bended rib- 
bon with reference to the ideal graphene sheet. We ob- 
served with good accuracy the validity of the relation 
= eis^ -I- e'fj , further proving that the total strain 
in a bended ribbon is the sum of two different contribu- 



tions: the first one (ssJ) is directly related to the bending 
process and the second one (els'*) is originated by edges 
effects, i.e. by the finite size of the nanoribbon. Although 
the first term seems to be quite negligible with respect 
to the second one, the previous energetic analysis reveals 
that both contributions are essential in order to explain 
the discrepancies between continuum and atomistic re- 
sults. 



IV. CONCLUSIONS 

In conclusion, we offered robust arguments suggest- 
ing that the correct value for the bending rigidity of a 
carbon nanoribbon corresponds to k = 1.40 eV, as cal- 
culated either through Ub = ^bended _ jrjfiat through 

Ub = E'^Xii, - i?f^;*,,,,,„,,. On the other hand, the 

relation Ut = E^^lf^H^ - E^'j^^^ leads to incorrect re- 

suits because of the emergence of a stretching field ■ 
We have further proved that such an in-plane strain field 
can be decomposed in a first contribution e^J'J due to the 
actual bending and a second one Ess due to the edges 
effects. 
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